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Abstract 

The production of D~ relative to Df as a function of xp with 600 GeV/c S~ beam 
is measured in the interval 0.15 < xp < 0.7 by the SELEX (E781) experiment at 
Fermilab. The integrated charge asymmetries with 600 GeV/c beam (0.53zb0.06) 
and 7r~ beam (0.06 it 0.11) are also compared. The results show the E~ beam 
fragments play a role in the production of D~, as suggested by the leading quark 
model. 
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1 Introduction 



Perturbative Quantum Chromodynamics calculations in leading or next-to- 
leading order predict very small or no asymmetry in the Xp or pt distribu- 
tions for charm and anticharm production[l,2]. However, fixed target data 
show some asymmetry between the production of some charm and anticharm 
hadrons in hadron-hadron interactions [3, 4]. SELEX has shown strong, beam- 
dependent asymmetries in production [5]. This experiment finds that Df 
production from a beam (but not a 7r~ beam) also has a large production 
asymmetry. This asymmetry could be due to the fact that the beam hadron 
shares a quark with one of the charge states (hence leading particle) and not 
with the other charge state (non-leading) . This is sometimes called "the lead- 
ing particle effect." 

For a ll^{sdd) beam the D^{cs) shares an s quark with the beam hadron and is 
a leading particle, whereas Df{cs) is not. For a n^ijid) beam, neither D~ nor 
is leading. Several theoretical models have been proposed to explain charm 
hadroproduction in the framework of non-perturbative hadronization. Among 
the proposed models are the color-drag string model[6], which is pronounced 
at high xp and is independent of pt, and the intrinsic charm model[7], which 
manifests itself at low pt and larger xp- 
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2 Appeiratus 



Data were taken during the 1996-97 fixed target run at Fermilab. The 600 
GeV/c negative hadron beam used in this measurement was composed of ap- 
proximately 50% E~ and 50% 7r~. Beam particles were tagged with a transition 
radiation detector system located upstream of the charm production target. 

A segmented target consisting of two copper sheets and three diamond sheets, 
each spaced by 1.5 cm, was used to produce charm particles. The total target 
thickness was 4.2% of an interaction length for protons. 

The SELEX experiment used a three-stage spectrometer designed for large 
acceptance at xp > 0.1 and for detecting the decay of charm particles. Each 
stage included a bending magnet and a detector system. 

SELEX used an online trigger to identify charm particles. The hardware trig- 
ger required at least 4 charged hadrons in the forward 150 mrad cone and 2 
hits from positive track candidates in a hodoscope after second magnet (M2). 
The software trigger made a full vertex reconstruction of the beam track and 
all tracks in the M2 spectrometer. Trigger conditions were also included in the 
simulation. 

The RICH detector, located after the second spectrometer system, was filled 
with neon gas at room temperature and 1.05 atm pressure. It identified charged 
hadrons whose trajectories went through the fiducial volume, typically requir- 
ing]? > 23 GeV/c[8]. Full reconstruction of the secondary vertex was provided 
by linking RICH-identified tracks back through the second stage magnetic 
spectrometer to the vertex silicon detector and associating them at a common 
decay vertex downstream of the primary interaction vertex. 



3 Data Analysis 

Initial data selection required two kaons to be identified by the RICH (ratio 
of the likehhoods: C{K)/C{'k) > 1) and the primary vertex to be in a target 
sheet. 

The criteria used to select Dg candidates included the following cuts: 

i) Secondary vertex separation significance L/a > 9 where L is the longitudinal 
separation between primary and secondary vertex and a is the error on L. 

ii) Secondary vertex was at least 100 microns outside of the target material. 
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iii) Each secondary track was extrapolated back to primary vertex ^-position 
to evaluate the transverse miss distance b. The second largest miss distance 
had to have {b/abY > 8, where ah is the error on b. 

These cuts were chosen to reject as many background events as possible with- 
out losing too much signal. They were optimized using real background and 
simulated signal events by maximizing the so-called significance: S/ \^Ns + 
where S was the yield from a Monte Carlo data set. The numbers of signal 
{Ng) and background (iVfe) events inside the square root were taken from data 
(all the events within the mass interval of 50 MeV/c^ centered at Dg mass 
value). The cuts are identical for the charge conjugate modes. None of the 
results presented here is sensitive to this optimization procedure. 

Since the RICH detector does not separate particles with absolute certainty, 
we expect some small amount of misidentification between pion and kaon 
that causes a refiection of D"^ under the Df peak (Detailed work on the 
contamination has been reported in Ref [9] measurement of Dg lifetime) . Only 
resonant ( 07r, K*K) channels were considered for this analysis to reduce the 
contribution of these reflections significantly. Df charm meson decays to 07r^ 
were selected by starting with candidate — >■ K~K^ decays. The invariant 
mass for two well-reconstructed oppositely charged tracks, identified as kaons 
by RICH, was calculated. candidates were those pairs whose invariant mass 
was within ±10 MeV/c^ of the mass (1020 MeV/c^). Similarly, those KKt: 
combinations that include a Xtt pair with an invariant mass value within 
±35 MeV/c^ of 7^°* mass (892 McV/c^) were selected as Df K^*{S92)K^ 
decays. We have obtained 172±14 D~ and 54±8 in 07r channel and 174±14 
D~ and 71±12 Df in K*K channel for xp >0.15 with the E~ beam. 

With the above cuts the Df peaks are clearly evident in the invariant mass 
spectra of KKt^~ and KKn'^ (Figure 1). This figure also shows Cabbibo- 
suppressed D"^ peaks. There is a clear excess of D~ over Df as seen in the 
figure. 



4 Production Asymmetry and xp Analysis 

Determining the yields through fitting the mass histograms for a specific xp 
value by a Monte-Carlo generated shape is often inaccurate for small statis- 
tics and fluctuating backgrounds. In SELEX the yield at a speciflc xp value 
was estimated using a sideband subtraction method assuming a linear back- 
ground. The mass ranges of the sideband background windows were [1.900 
GeV/c2, 1.940 GeV/c^] and [2.060 GeV/c^, 2.140 GeV/c^]. We defined asym- 
metric sidebands to avoid the infiuence of D"^ —>■ K^K~t{^ and to exclude 
the D*(2010) mass region. 
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Fig. 1. Invariant mass distributions of (a) KKtt^ and (b) KKtt^ from the S~ 
beam. Mass interval and sidebands used in determining the yields for asymmetry 
calculations are shaded in both histograms. The data shown here include events 
having the decays — > (f>TT^ or 



The yields after sideband subtraction were corrected for the acceptance (recon- 
struction efficiency and geometrical acceptance) of the detector. To estimate 
the acceptance, Ds events were generated by a Monte Carlo program with a 
flat distribution in and a Gaussian-distributed transverse momentum with 
mean pt=0.8 GeV/c. In a given simulation data set the Dg decays only into 
the K*K or (pn mode. Decay tracks were digitized after smearing with de- 
tector resolution and multiple Coulomb scattering effects. The detector hits 
were OR'ed into the hit banks of interaction data. The new hit banks were 
passed through the SELEX off-line software. The acceptance was measured 
using the ratio of the reconstructed events to the embedded events. The set 
of cuts that was used to extract the signal was applied in this case as well. 
The most important issue for the asymmetry measurement was the relative 
efficiency of Dj and Df. 



As indicated in Figure 2, the average difference in acceptance between D~ and 
is very small over all xp range compared to the statistical uncertainty. 
This shows that the spectrometer is charge independent for decay events. 
In figure only the acceptance for (pn channel is shown. The acceptance for 
K*K channel is slightly lower. The acceptance in xp is independent of pt- 



5 Results 



After all the cuts and acceptance corrections, the resulting Dg yields as func- 
tions of xp are shown in Figure 3 and listed in Table 1. Resulting data points 
are fit to a functional form (1 — xp)^- The values of n obtained from the fits 
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Fig. 2. Acceptances for D~ and Df obtained by embedding Monte Carlo events 
into data events. The acceptance here is the combination of geometrical acceptance 
and reconstruction efficiency of the ^tt channel. 

are shown on the figure. 
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Fig. 3. Acceptance-corrected xp distributions for D~ and Df from beam. Fits 
of the yields to (1 — xpY^ for each charge state arc plotted and the n- values listed. 

Since the beam flux was the same for Df, the two distributions compare the 
relative production differential cross-sections for these states. Figure 3 shows 
that, for the E~ beam used in this measurement, production is favored 
over Df at all xp, and the difference increases at large xp. 

We can discuss this difference in terms of an asymmetry parameter A, defined 
as 



where Nj^- and Njj+ are the corrected yields for and D^, respectively. The 
asymmetry was calculated for five equally divided bins over a.n xp range of 
0.15 to 0.40 and for three equally divided bins over an xp range of 0.40 to 0.70. 
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XF 






Asymmetry 


0.15-0.20 


875±156 


619±141 


0.17±0.14 


0.20-0.25 


669±102 


409±91 


0.24±0.13 


0.25-0.30 


723±92 


262±63 


0.47±0.11 


0.30-0.35 


413±62 


125±43 


0.54±0.13 


0.35-0.40 


282±49 


69±38 


0.61±0.18 


0.40-0.50 


232±34 


47±26 


0.66±0.16 


0.50-0.60 


124±28 


5±7 


0.92±0.10 


0.60-0.70 


15±9 


5±5 


0.50±0.43 



Table 1 

Summary of D~ and Df yields and asymmetries from S~ beam. The errors are 
statistical only. Yields are obtained from resonant state K*K and 07r events. 

Figure 4 displays the acceptance-corrected asymmetry as a function oi xf for 
the S~ beam. It shows that there is a significant asymmetry in favor of D~ 
and that the asymmetry increases gradually as increases. The asymmetry 
values are included in Table 1. 



I 0.9 r 
E 0.8 r 




^ — ' — • — • — • — ' — • — • — • — ' — • — • — • — ' — • — • — • — ' — • — • — • — I 

0.2 0.3 0.4 0.5 0.6 0.7 

Fig. 4. Production asymmetry for beam data as a function oixp- Yields obtained 
from resonant {K*K and (j)!^) events were used to calculate the asymmetry. 



Figure 5 displays the asymmetry as a function of p^. These asymmetry values 
are also corrected for acceptance. One can see that the asymmetry is flat 
within the observed range (up to pi < 5GeV/ c) 

In order to further explore the leading particle effect, the analysis was repeated 
with data from the vr^ beam, obtained under the same conditions as with the 
S~ beam. Because pion beam statistics are much lower, we compare only 
integrated asymmetry results for all KKn events with xp > 0.15. For the 7r~ 
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Fig. 5. Production asymmetry for S~ beam data as a function ofp^. Yields obtained 
from resonant {K*K and (j)ir) events were used to calculate the asymmetry. 

beam, which has no leading particle, the integrated asymmetry is consistent 
with zero (A = 0.06±0.11). On the other hand, analysis of the beam data 
in the same way results in a large asymmetry in favor of D~ {A = 0.53±0.06). 



6 Systematic Study 

Studies of possible systematic errors due to the side-band subtraction method 
were done by varying the size and position of the side bands. Effects of chang- 
ing xp bin sizes on the results were also investigated. Systematics of the ac- 
ceptance calculations have been checked with meson asymmetries and polar- 
izations, all of which should be zero and are. The false asymmetry due to 
our hardware trigger was also studied. Even when combined together, these 
effects are all considerably smaller than the statistical uncertainty and are 
neglected. The contribution of the misidentification of tt" beam particles as 
E~ is estimated to be a few percent. The resulting dilution in the asymmetry 
is negligible. 

As mentioned before, background including the D"^ contamination under the 
Df peak is highly reduced by limiting data to the resonant states. Effects of the 
remaining background were studied by comparing the integrated asymmetries 
obtained from the two resonant states. In the 07r case all backgrounds, includ- 
ing contamination, are negligible and cannot affect the asymmetry. The 
07r integrated asymmetry is 0.52±0.06. For the K*K channel it is 0.42±0.08. 
The effect of the contamination reported in ref [9] would reduce a real K*K 
asymmetry of 0.52 to an observed value of 0.48, statistically compatible with 
observation. The overall dilution effect is much smaller than the statistical 
uncertainties in individual bins and is not included in the final results (Table 
I and Figure 4). 
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7 Conclusion 



To summarize: the 'E~{dds) beam data show a strong production asymmetry 
favoring Dj{cs) production. This is consistent with leading particle effects. 
However, the integrated asymmetry from 7r~ beam at xp >0.15 for Df meson 
is 0.06±0.11, which is consistent with zero asymmetry as expected since nei- 
ther nor D~ is a leading particle. Our E~ results are consistent with the 
previous measurement done by WA89 experiment at CERN with 340 GeV/c 
E~ beam. 

The SELEX pion result of negligible integrated asymmetry agrees with the 
higher-statistics differential distribution for a 500 GeV tt" beam reported by 
E791[4]. Their integrated asymmetry in the xp range 0.1 to 0.5 is 0.032±0.022. 
Our results also favor the color drag model over the intrinsic charm model, 
since the color drag model predicts a large asymmetry at large xp independent 
oipt[A,10,ll]. 
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